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Abstract 

We examine the feasibility of deriving neutrino mixing parameters 6 and ^13 from 
the cosmic neutrino flavor composition under the assumption that the flavor ratios of 
the cosmic neutrinos at the source were F^^ + Fp^ : F^^ + Fp^ : F^^ + Fp^ = 1:2:0. 
We analyze various uncertainties that enter the derivation of 6 and 613 from the ratio 
of the shower-like to /i-tracks events which is the only realistic source of information 
on the flavor composition at neutrino telescopes such as ICECUBE. We then examine 
to what extent the deviation of the initial flavor ratio from 1:2:0 can be tested by 
measurement of this ratio at neutrino telescopes taking into account various sources of 
uncertainty. 
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1 Introduction 



According to the current models, the astrophysical objects such as sources of Gamma Ray 
Bursts (GRBs) [1], type I b/c supernovae pj and Active Galactic Nuclei (AGNs) [3] can 
emit beams of neutrinos luminous enough to be detectable at the neutrino telescopes that 
are under construction. The AMANDA experiment [1] at the south pole, which came to the 
end of its mission in 2006, has set the following bound on the diffuse flux of neutrinos 

< 8.2 GeV cm-^ sr"^ yr'^ . (1) 

A km^-scale neutrino telescope named ICECUBE is under construction which encompasses 
AMANDA. If the bound in (P is saturated, the completed ICECUBE can collect ~ 4000 
cosmic neutrino signal each year [5]. Notice that this bound is on the sum of the neutrino 
fluxes from sources at cosmological distances. In principle, core collapse supernova explosions 
leading to an intense neutrino flux detectable at km^-scale neutrino telescopes can also take 
place in the close-by galaxies located at a distance of < 10 Mpc [HIE!- If such an explosion is 
registered during the time that the ICECUBE is in full swing, ICECUBE can record about 
a few hundred neutrino events from a single explosion [7j. In addition to ICECUBE in the 
south pole, three neutrino telescopes NEMO [8], ANTARES [9] and NESTOR [lO] are being 
constructed in the Mediterranean sea. Moreover, the so-called KM3NET neutrino telescope 
[TT] is planned to be constructed in the Mediterranean sea. 

In view of this prospect, extensive studies have been performed on the possibility of 
deriving information on the mixing parameters of neutrinos by studying the flavor ratio of 
the neutrinos at the detector [121 [131 [HI- The method is based on the following argument. 
Suppose the flavor ratio at the source was F^^ + Fp^ : F^^ + Fp^ : F^^ + Fp^ = : w^^ : . 
After propagating the distance between the source and the detector, the flavor ratio will 
become 

F.e + F,^ ■■ F., + Fp, : F,^ +F,^ = J2 <Po^e : J] ^P^^ : J] wIP^^ , (2) 

a a. a 

where Pq,/3 is the probability of — * vp- Considering the very long distance between the 
source and the Earth [i.e., Am^jL/ {2E^) ^ 1), the oscillatory terms in P^/? average out ^: 

= P(z/„ ^ up) = P(z7, ^up) = Y^ (3) 



•^For a detailed discussion of the loss of coherence, see 
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where Uai are the elements of the neutrino mixing matrix. Notice that Pap is independent 
of the neutrino energy E^, distance L and the mass square differences Am^^. 

In a wide range of models that lead to detectable cosmic neutrino flux, the neutrino 
production takes place through t:^ — > ^ and subsequently /i e^Vf.v^. The flavor 

ratios at the source are therefore predicted to be : w° : w° = 1 : 2 : 0. Thus, by 
measuring the flavor ratio at Earth, one can in principle derive the absolute values of the 
mixing matrix elements which yield information on the yet-unknown neutrino parameters 
6*13 and 5 as well as the deviation of 6*23 from 7r/4 [121 [131 El IE]- It is also suggested to 
employ cosmic neutrinos to discriminate between the standard oscillation scenario and more 
exotic possibilities [TTl [18] . 

The flavor identification power of ICECUBE and other neutrino telescopes is limited. 
In fact, in the energy range of interest (100 GeV < < 100 TeV), ICECUBE can only 
distinguish between shower-like events and the /i-track events. Each of these two types of 
events can receive contributions from different flavors. As discussed in the next section, 
several input and assumptions go into derivation of the flavor composition from the ratio of 
the shower-like events to the //-track events. Lack of knowledge or uncertainty on these input 
parameters will lead to uncertainty in derivation of 613 and 6 from the cosmic neutrinos. To 
derive information on the unknown mixing parameters (^13 and/or 6) from the ratio, one 
should also consider the effect of the uncertainty on the known neutrino mixing parameters, 
623 and 6i2- A complete treatment of all these effects is missing in the literature. In the first 
part of this paper, we study the possibility of deriving ^13 and S from the ratio of shower-like 
events to /i-track events. We take into account the uncertainty of the aforementioned inputs 
as well as the possible uncertainty in the measurement of the ratio itself. In our analysis, we 
include contributions to the shower-like and /t-track events whose effects can be larger than 
the effect of 6 and ^13 but their effects have been overlooked in the literature. We discuss 
how much the precision of various input parameters has to be improved in order to make 
the measurement of 613 or 6 by neutrino telescopes a reality. 

Various effects can cause a substantial deviation of the flavor ratio from : : = 
1:2:0. In the second part of the paper, considering the realistic uncertainties in the input, 
we discuss to what extent the ratio at the source can be determined under the assumption 
that propagation of neutrinos from the source to detector is simply governed by the standard 
oscillation formula in Eqs. (!2]l3l) . 

The paper is organized as follows. In Sect. [21 the mechanism for flavor identification 
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at neutrino telescopes is described. In Sect. [3l the features of the cosmic neutrino fluxes 
predicted by the mainstream models are discussed. In Sect. HI assuming CP conservation, 
the effects of possible sources of uncertainties on the derivation of S13 are studied. An analysis 
of derivation of 6 in the presence of uncertainties is performed in Sect. [51 Sect. [6] is devoted 
to discussing various possible mechanisms through which the initial ratios can deviate from 
1:2:0. We discuss whether by measuring the flavor ratio on Earth, it will be possible to 
differentiate between models. A summary of the conclusions is provided in Sect. [3 



2 Flavor identification 

ICECUBE and its Mediterranean counterparts can basically distinguish only two types of 
events: 1) shower-like events; 2) yU-track events. As shown in the seminal work by Beacom 
et al. [20] , one can derive information on the flavor composition by studying the ratio of the 
yU-track events to shower-like events. Let us define 

^ Number of Muon-track events 
Number of Shower-like events 

There is a threshold energy, Eth below which the neutrino cannot be detected by a neutrino 
telescope. The value of Eth depends on the structure of the detector and the type of the 
event (shower-like versus /i-track). Since the detection of neutrinos coming from above suffers 
from the large background from cosmic rays, the neutrino telescopes mainly focus on upward 
going neutrinos; i.e., neutrinos that pass through the Earth before reaching the detector (see 
however, [2T]). The Earth is opaque for neutrinos with energies higher than ~ 100 TeV [22] . 
In sum, the neutrino telescopes mainly study the muon neutrinos in the range Et^ ~ 100 GeV 
and E*^^ ~ 100 TeV. In calculating the fluxes we set the upper limit of the integration equal 
to El^* = 100 TeV. By using this value for the upper limit, we can neglect the attenuation 
of the neutrino flux crossing the Earth which depends on the direction of the neutrino [23] . 

Two sources contribute to the /x-track events: (i) Charged Current (CC) interaction of u^^ 
or z/^ producing /i or ft; (ii) CC interaction of Ur and Ur producing r or r and the subsequent 
decay of r and f into fi and /2. In the literature, the contribution of z/^ (via t fi) 

to /i-track events has been overlooked but to study the effect of 613, one should take into 
account such sub-dominant effects. 
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The contribution of and z/^ to the /i-track events can be estimated as 



pANa II R,iE,,E^,, 



dFy„ da 



cc 



dEj,^ dEfj_ 



■dEfj^dEu^ + [particle — > antiparticle], 



(5) 



where p, A and A^^ are respectively the density of the medium (ice/water), the effective 
area of the detector and the Avogadro number. R^{Ei, E2) is the muon range which is the 
distance traveled in the medium by a muon with energy Ei before its energy drops below 
E2. The muon range in ice is given by [24| 



R^{E,,E2) = (2.6 Km) In 



2 + 4.2 X 10-3^1 
2 + 4.2 X 10-3^2 



(6) 



where both Ei and E2 are in GeV. Finally, dF^^j dE^^ and a^'" are respectively the neutrino 
flux spectral function at the detector and the cross section of the CC interactions of z/^. 

The contributions from Vr and Vr to /i-track events can be estimated as 



BpANA 




dF^^ da^^ 
dEy^ dEj. 



f {E^, E^) R^{E^, E^^)dE^dE^dE,^ + (z/, ^ z/,) 



(7) 



where B = Br(r pv^Vr) = 17.8%. The function f{Er,Efj^) in the above equation is the 
probability density of the production of a muon with energy E^ in the decay of a r lepton 
with energy Er- That is 



f{Er, E^ 



1 dViriEr) piE^)u^u^) 



r 



dE„ 



The details of the calculation of fiE^-, E^) can be found in Appendix Rl 

Three types of events appear as shower: i) the Neutral Current (NC) interactions of all 
kinds of neutrinos; ii) the CC interactions of and Ve'-, hi) the CC interactions of Vj. (z/,-) and 
the subsequent hadronic decay of r (f). Showers from NC interaction of all three neutrino 
flavors can be estimated as 



l=e,ii,T 



(9) 



where L is the length of the detector. The rate of the electromagnetic showers from the CC 
interactions of z/p and z/p is 



pALN/ 



dF, 



^cc 



dE,, 



a'^'-dE,^ + 



dF, 



v„ -cc 



dE„ 



(10) 
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The rate of showers originated from the CC interaction of Ur with the subsequent hadronic 
decay of the r lepton is 



'1-B)pALNa 



Ecut J T? fEcut J Tf<_ 



Notice that while the shower-hke events are given by the length of the detector, the /x- 
track events are given by the muon range [see Eqs. (15|71) ]. This is because muons can emit 
Cherenkov hght and trigger the detector even if they are produced outside the detector but 
within the range i?^. In other words, for the muon detection, the effective volume is larger 
than the geometrical volume. 

To write down the above formulas, several simplifications have been made: 

• Obviously, neutrinos entering the detector through different zenith angles have prop- 
agated different lengths inside the earth so the amount of attenuation is different for 
them. In other words, to be precise, the zenith angle dependence of E'^^ has to be 
taken into account. 

• The energy threshold for detecting the neutrino also depends on the direction. For the 
vertically propagating muon at ICECUBE, E^^^ is about 20 GeV; that is while, for the 
muons propagating horizontally, E^^ is about 100 GeV [25j . 

• For the high energy muons with E^> 1 TeV, the muon ranges both in the ice and rock 
exceed 1 km, (see Eq. © and Ref. [21]). The depth of the ice at the site of ICECUBE 
is about 2810 m; so a considerable number of muons reaching the ICECUBE would be 
produced in the rock beneath the ice where the density is quite different. 

• As mentioned before, for neutrinos with energies higher than 100 TeV the Earth is 
opaque. For very high energies, Vr can however be regenerated through z/,- — > r ^ 
Vt- ... Vr. As a result, and with Ey ^ 100 TeV can give a contribution to 
the upward-going neutrino flux with E < 100 TeV. The neutrino flux at high energies 
is expected to be suppressed. Thus, such a contribution is expected to be negligible 
|26] . This assumption can in principle be tested by measuring the downward-going 
shower-like events (which have not traversed the Earth). 

Throughout the present analysis, we use the approximate formulae (15]I71I9][T0|TT|) . We deter- 
mine how much the uncertainties in various inputs have to be improved in order not to be an 
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obstacle for determination of S and S13. As we shall see, our conclusion is that even without 
the above subtleties, the required precision in certain input parameters is so fine that seems 
beyond reach in the foreseeable future. Taking into account the above uncertainties not only 
will not change our conclusion but will further confirm it. 



3 The standard picture 

To evaluate R, several input parameters have to be known: i) the energy spectrum of the 
incoming neutrinos; ii) the ratio of the neutrino flux to the anti-neutrino flux; iii) the initial 
ratio : w° : w^. We rely on the predictions of the models for such input parameters. 
Although the models differ in details, they share some common features. From now on, 
we call these features the "standard picture." The features of the standard picture are 
enumerated below. 

• In the standard picture neutrinos are produced in the following chain of processes. 
First, the energetic protons in jets collide on 7 or on the background protons and 
produce vr"*" and vr". Then, 



7r+ /i+Z/^ TT ^ jj, i^^ 



(12) 



Thus, w° : w° : w° = 1 : 2 : 0. 

The energy spectra of the neutrinos follow power law distributions: 

where Mp is the normalization factor for each neutrino and anti-neutrino flavor, a is 
the spectral index. In the standard picture where the initial protons are accelerated to 
high energies via Fermi acceleration mechanism, the spectral index is expected to be 
equal to 2 [27] . 

Regardless of the relative amount of 7r+ and tt^, we expect = Afp^. However, 
AfpJAfu^ depends on the initial composition of 7r+ to which in turn depends on the 
details of the model. 
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4 Uncertainties and their impact on ^13 measurement 



Our knowledge of the sources of the cosmic neutrinos is quite hmited and mostly speculative. 
A myriad of known and un-known effects can cause deviation of the initial flux from the 
standard picture that was described in the previous section. In this section, we compare 
the effect of a deviation from the standard picture on R with the effect of a nonzero S13. 
Here, we assume that the neutrino mass matrix conserves CP. A discussion of CP- violation 
is given in sect. [H 

Let us fix our convention for the mixing angles. Here, we use the standard parametrization 
of PDG [2H] for the neutrino mixing matrix with [29] 

< 013 < 0.2 < ^ and < 5 < 27r . 

The sensitivity of Pa/s on the phase 6 is through sin 6*13 cos 5; so the CP-conserving cases with 
6 = and 6 = n will have distinct effects. We will consider both cases. Throughout this 
section we take : w° : w° = 1 : 2 : 0. 

In sect. 14. 1^ we discuss the effect of the variation of the spectrum on ratio R. In sect. 14. 2^ 
we discuss the dependence of R on Afy^/Nu^. Sect. 14.31 gives a brief discussion of the neutrino 
nucleon uncertainty and its effects. 

4.1 Energy spectrum of incoming flux 

As mentioned before, in the standard picture, the neutrino flux follows a power-law spectrum 
of form Eq. f[T^ with a = 2. However, more careful considerations of the details of the Fermi 
acceleration and the properties of the target particles show that a can deviate from 2 and 
take any value in the range (1,3) [161 l30l [3T] . 

The energies of the muons and showers entering a neutrino telescope can be measured. 
However, extracting the energy of the incoming neutrinos that induce such events is not 
straightforward. In the case of /i-track events, the muon can lose a substantial part of its 
energy before entering the detector. On the other hand, limiting the analysis to the muons 
produced inside the detector will reduce the statistics. In the case of the shower-like events 
originating from the NC interaction of neutrinos, the energy of the shower does not give the 
energy of the initial neutrino because a part of the energy is carried away by the final neutrino 
which escapes detection. Nevertheless, it is shown in [2Qj that for Ef,dFy/ dE,j = 0.25 GeV 
cm~^ sr~^ yr^^ after one year of data-taking, a can be determined with 10 % uncertainty. 
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Figure 1: The dependence of R on sin^ ^13 for different values of the spectral index, a. The 
thicker lines correspond to 5 = tt and the thinner ones correspond to 5 = 0. We have used 
the central values for the neutrino-nucleon cross section [22] and have set Np^/Ny^ = 0.5 and 
(A/p^ + ■N'ufj) / {J^ue + -^ue) — 2- The input for 6*12 and 6*23 are set equal to the best fit in [29j. The 
vertical line at 0.041 shows the present bound at 3(7 [29]. 

Fig. [1] shows R versus sin^^is for cos 5 = ±1 and various values of a. As seen from the 
figure when 5 = 0, the sensitivity of R to 5^3 is very mild and less than 2 %. That is while 
for cos 5 = —1, the sensitivity to s^g is about 10 %. The disparity between cos 5 = +1 and 
cos 5 = —1 means that for 5^3 ~ 0.04, the contributions from S13 cos 5 and are comparable. 
In fact, expanding R in powers of S13 confirms this claim: 

i? ~ ri + r2Si3 cos(5 + r3S^3 cos^ 5, (14) 

where for the central curve with a = 2.0, ri ~ 2.55, r2 — —0.66 and ~ 2.65. If 6'23 deviates 
from 7r/4, in addition to the sl^cos^6 term, a term proportional to sf^ has to be added to 
Eq. (m. 

As seen from the figure, even for cos 5 = —1, the sensitivity to 5^3 can be obscured by 
the 10 % uncertainty in a. However, for 5^3 > 0.02, the bands between a = 2.2 and 1.8 for 
cos 5 = 1 and cos 5 = —1 have no overlap. This means that for sf^ > 0.02, 10 % precision in 
a is enough to distinguish cos 5 = 1 from cos 5 = — 1. 
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Notice that the curve with a = 2 is closer to that with ot = 2.2 than that with a = 1.8. 
This means that the effect of the uncertainty decreases by increasing the value of a. 



R 




Sin^e 
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Figure 2: The dependence of R on sin^ ^13 for different values of the p parameter defined in Eq. (llSp . 
We have set 6 = n. The rest of the input parameters are the same as in Fig. [TJ The vertical line 
at 0.041 shows the present bound at 3a [29]. 



Considering the unknown nature of the production mechanism, it is not dismissed that 
the energy spectrum of neutrinos does not follow a simple power-law form. For example 
the spectrum can be a sum of two power-law functions each originating from a separate 
mechanism {i.e., aE~°'^ + bE~°'^). We study such a possibility in Fig. [21 where we have taken 
the shape of the spectrum to be of the form: 



E 



-1 



100 TeV 



(15) 



where p is a dimensionless parameter that determines the magnitude of the second term. 
Both terms can originate from the Fermi acceleration mechanism [32] . Notice that the second 
term is subdominant. The curves from up to down respectively correspond to p = 1, 0.5, 
0.3, 0.1 and 0. As seen from the figure, the uncertainty on p obscures the extraction of the 
mixing angle 613. 
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4.2 Uncertainty in HyJHy^ 

Since the source is made of matter rather than anti-matter, we in general expect tt^ to 
dominate over vr" and therefore < A/p^/A^^ < 1. There is not any estabhshed or pro- 
posed method for determining M^^jMy^ in the neutrino telescopes in the energy interval 
(100 GeV, 100 TeV). As a result, this ratio appears as a source of uncertainty in determina- 
tion of R. Curves in the Fig. [3] show the dependence of R on sin^ 6*13 for two extreme values 
A = My^jMy^ = and A = 1. For 6 = n, the variation of Afp^/Af^^ in the interval [0, 1] causes 
a change in R of about 5 % which can obscure the determination of S13. Notice that for 
> 0.005 the bands between A = 1 and A = for 6 = and 6 = it are separate, so the 
uncertainty in A will not cause a problem for discriminating between cos 5 = ±1. 



R 
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Figure 3: The dependence of R on sin^ ^13 for different values of the parameters A = Mv^jMy^. 
The thicker lines correspond to (5 = vr and the thinner ones correspond to (5 = 0. The spectral index 
has been set equal to 2. The rest of the input parameters are the same as in Fig. ([T|. 



4.3 Uncertainties in cross sections 

To calculate the cross section, information on the Parton Distribution Functions (PDFs) of 
the nucleon is needed. The center of mass energy of a system composed of a neutrino with 
energy ~ 100 TeV incident on a proton at rest is {^E^m^^l'^ = 450 GeV. The center 
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of mass energy of the e — p HERA collider is about 320 GeV. Thus, to calculate a^N in 
the energy range relevant for this study {i.e., 100 GeV < 100 TeV) the results of the 

HERA experiment can be employed. The current uncertainty on the PDFs is about 3 % 
[33] . LHC can further improve the precision of the PDFs. 

The uncertainty in the cross section of all types of neutrinos (each flavor of neutrino and 
anti- neutrino) originates from the same uncertainties in the PDFs. As a result, the resultant 
uncertainty in the numerator and denominator of R cancel each other so the certainty in the 
cross section will not be a limiting factor for determining ^13 (or 5) from the cosmic neutrino 
flavor composition. 

5 Determination of 5 

In the literature it is suggested to derive the Dirac CP-violating phase {5) by studying the 
flavor composition of the cosmic neutrinos [T3| [Hj. Considering the expenses and chal- 
lenges before measuring this phase through the more conventional proposals {i.e., neutrino 
factory or superbeam methods), it is worth giving this possibility a thorough consideration. 
However, in the literature the flavor identiflcation power of the detector has not been real- 
istically treated. To be speciflc, the quantities that have been previously analyzed in the 
context of deriving 5 are ratios such as R' = Fy^/{Fy^ + F^^) [T2[ [13], E] which cannot be 
directly derived at neutrino telescopes. In this section, we assess the possibility of measuring 
5 considering realistic flavor identiflcation power of neutrino telescopes {i.e., studying ratio 
R) and taking into account various sources of uncertainty for the flrst time. 

By the time a statistically signiflcant number of cosmic neutrino events is collected, we 
expect noticeable improvement in determination of the input parameters. In particular, in 
the case that the parameters are in favorable range, we expect progress in the following 
measurements: 

• If sin^ is close to the present bound, the forthcoming experiments can measure its 
value with a precision of Asin^ 6*13: 

sin^ = sin^ ^13(1 ± Asin^ 6^/ sin^ ^13). 

In fact, for relatively large values of sin^ 61^, the uncertainty Asin^ 61^ j sin^ ^13 can be 
reduced to as small as 5 % [311 ES] . 
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• In case that statistically significant number of cosmic neutrinos are collected, R can 
be measured with an uncertainty of AR: 

R = R{1±AR/R). 

As shown in [20|, a precision of AR/R ~ 7% can be obtained provided that the 
number of events exceeds ~ 300. This would be achieved with a neutrino flux of 
E^dF^/dE^ = 0.25 GeV cm~^ sr^^ yr^^ after a couple of years of data-taking. 

• The forthcoming long-baseline [36] and reactor neutrino [37] experiments can respec- 
tively measure the solar and atmospheric mixing angles by a precision of ~ 6 %. That 
is 

sin2^i2 = sin2^i2(l±6%), 
sin^^ss = sin2^23(l±6%). 
The present best-fit values are sin^ 9i2 = 0.32 and sin^ 623 = 0.5. 

Remember that the ratio Afp^/Afu^ cannot be measured. Since the initial jets creating 
the charged pions (and subsequently the neutrinos) are mainly made of protons rather than 
anti-protons, we expect My^jMy^ < 1. Considering all these uncertainties, the question is 
whether it will be possible to extract the value of 6. 

Fig. m addresses this question. Drawing the plot, we have assumed that R will be found to 
have a typical value of 2.53 with an uncertainty of AR/R. This value of R can be obtained 
by taking maximal CP-violation {6 = vr/2), sin^ 6^13 = 0.03, : w° : w° = 1 : 2 : 0, 
a = 2 and My^/Mv^ = 0.5 . We have looked for solutions in the 5 — a plane for which 
R = 2.53(1 ± AR/R), varying the rest of the relevant parameters in the ranges indicated in 
the caption of Fig. |H The regions covered with dots, little triangles and crosses respectively 
correspond to 7%, 1.5% and 1% precision in the measurement of R. Notice that the figure is 
symmetric under 6 27t — 6. The symmetry originates from the fact that the dependence 
of i? on (5 is through cos 5. As mentioned earlier, a can be independently determined by 
the measurement of the energy spectrum with about 10% precision. (For a = 2, the direct 
measurement of the energy spectrum can restrict the value of a to the region between the 
vertical lines at a = 1.8 and 2.2.) As seen from the figure, with AR/R = 7%, 6 cannot 
be constrained. In fact, any point between the vertical lines can be a solution. The figure 
shows that reducing AR/R to 1 % (but keeping the rest of the uncertainties as before), some 
parts of the solutions can be excluded. In particular, the region around 6 = it will not be a 
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1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

a 

Figure 4: Points in the (a, S) space consistent with R = 2.53ibAii. True values of the (a, 6) pair are 
(2,7r/2). Points displayed by dots, triangles and crosses respectively correspond to AR/R = 7%, 
AR/R = 1.5% and AR/R = 1%. To draw this figure we have varied sin^ 6^13 e (0.028,0.032), 
sin2 0i2 G (0.30,0.34), sin^ ^23 G (0.47,0.53) andMpJAfu, G (0,1). 

solution anymore. Notice that along with 5 = 7r/2, 5 = is also a solution. This means that 
despite maximal CP-violation, the CP-violation cannot still be established. We examined 
the robustness of this result. We found that reducing the uncertainties in the mixing angles 
(even in ^13) will not noticeably change the overall conclusion. However, the sensitivity to 
AR/R seems to be high. Notice that with a precision of AR/R = 1.5%, there are some 
regions of solutions (covered by the triangles) that can be excluded if the uncertainty is 
reduced to 1% {i.e., they are not covered with crosses). As we will see below, the sensitivity 
to power index is also high. 

Drawing Figs. ([5]-aj5]-b), we have respectively taken R = 3.3 (corresponding to a = 1) 
and R = 3.59 (corresponding to a = 3). The rest of the input is the same as Fig. (j4]). 
From Fig. (l5]-b) we observe that for a = 3, the measurement of R with 7% uncertainty 
determines a with better than 6 % precision which will probably be more accurate than the 
direct determination of a from the energy spectrum measurement. Notice that when a = 1 
or 3, even a precision of AR/R = 1% will not be enough to constrain 6. However, from 
Fig. (l5]-a) we observe that in case of a = 1 and 6 = tt/2, if the error in direct measurement 
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Figure 5: The same as Fig. H] except that in figure (a), we have taken R = 3.30 (corresponding to 
a = 1) and in figure (b), we have taken R = 3.59 (corresponding to a = 3). Points displayed by 
dots and triangles respectively correspond to AR/R = 7% and AR/R = 1%. 

of a is reduced to 10% {i.e., if a is constrained to the region between the dashed vertical 
lines) and if R is measured with 1% precision, one can exclude solutions around 6 = n but 
CP-violation cannot still be established. 

As emphasized in [13j, the sensitivity of R to the variation of 623 is significant. In fact, the 
present uncertainty (i.e., sin^ 6*23 = CSOlg^g at 3a c.l. |29]) leads to a sizeable uncertainty 
of ~ 20 % in R. As mentioned earlier the forthcoming experiments [37] can improve the 
precision of sin^ 6*23 to 6 %. A variation of 6% in sin^ ^23 leads to a ~ 4 % change in R 
which is comparable to the effect of cos 5 for sin^^ia = 0.03. To pinpoint the effect of the 
uncertainty of 6'23 on R, we have presented Fig. [61 The input of Fig. [6] is similar to Fig. H] 
but to study the effect of 6*23, we have fixed AfpJAf^^, 612 and 613 to their central values. 
Comparing Figs. [Hand Ej we find that improving the precision of sin^ 6'23 from 6% to 1% can 
help us to remove a substantial part of the spurious solutions. Especially for AR/R = 1% 
solutions with Stc/A < 6 < 57r/4 can be removed by reducing the uncertainty of sin^ 6'23 to 
1 %. 
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Figure 6: Points in the (a, 6) space consistent with R = 2.53 it AR. As in Fig. [U the true values 
of the (a, (5) pair are (2,7r/2). Points displayed by dots and triangles respectively correspond to 
AR/R = 7% and AR/R = 1%. To draw this figure we have fixed sin^ (9i3 = 0.0301, sin^ 6*12 = 0.32, 
MpJMu, = 0.5 and varied sin^ 6^23 in 0.5(1 ± 1%). 

6 Initial flavor composition 

In the previous sections, in accordance with the "standard picture", we had assumed that 
: w° : w° = 1 : 2 : 0. Various mechanisms can intervene to cause a deviation from this 
simple prediction. For example, in the so-called stopped muon scenario, the muon production 
takes place inside a high magnetic field so the muons come to rest before decay. As a result, 
the neutrinos produced from the decay of the muon would be below the detection energy 
threshold. This effectively leads to : w° : w° = : 1 : [30]. On the other hand, a 
contribution from neutron decay, n i^eP^~, would increase w^/w'^. Even in the standard 
picture, along with the production of vr^, charged Kaons can also be produced in the jets. 
Like the case of the charged pions, the main decay mode of charged Kaon is ^ fiu^, so 
neglecting the subdominant modes Bt{K — > n^^v^) = 3.3% and Bt{K — > vr^ei^e) = 5% [28j, 
we would naively expect that flavor composition of the neutrinos from the Kaon chain also 
follow the standard picture (i.e., : w° : w° ^ 1 : 2 : 0) but there is a subtlety here. In the 
chain process, tt fiu^, /i — > z/gZ/^e, the average energy of each of the produced neutrinos in 
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the rest frame of the pion is about m,r/4. That is while in the case of Kaon decay, the average 
energies of the neutrinos are different: The energy of z/^ produced directly from the Kaon 
decay in the Kaon rest frame is {mj^ — m^)/(2mi<-) ~ 236 MeV; that is while, the average 
energies of the neutrinos from the secondary muon are [{mj^ + mp/(2mi^)]/3 ~ 86 MeV. 
As a result, limiting the detection to an energy range (e.g., (100 GeV,100TeV)) the ratio 
: ^0 {i.e., j^'^^\dFl^JdE + dFl/dE)dE : j^'^^\dF^JdE + dF^JdE)dE) will deviate 
from 1:2. In fact, the ratio : w° would depend on the energy spectrum of the initial 
charged Kaon. 

The mechanisms that we pointed out above are all processes that expected to exist and 
play at least a subdominant role within the framework of the mainstream models. None 
of these mechanisms creates r neutrino or anti-neutrino at the sources {i.e., they all yield 

: 1/7° : = : w° : 0). However, more exotic mechanisms can be at work to create v-r or 
Vr at source: In principle, the collision of pp or p^ at the source can create D meson whose 
decay can produce z/^ and Vt- According to [38j, the contribution of the D meson to the 
neutrino flux becomes important only for Ejj > 10^ TeV. For lower values of energy in which 
we are interested in, the contribution is about three orders of magnitude below the present 
bound. To our best knowledge, within the Standard Model (SM) of particles, in the energy 
range 100 GeV-100 TeV, remains much smaller than and w^. However, if neutrinos 
have some yet unexplored properties beyond SM, tau neutrinos or anti-neutrinos can be 
produced at the source. For example, suppose neutrinos posses tiny transition magnetic 
moments of form 

HreFapT^lCa'^'^i'r aud/or fi^^F^pu'l^Ca^'^iyr . (16) 

There are strong bounds on the transition moments, l/ia/j], from different terrestrial experi- 
ments and astrophysical observations [28l [39] : 

/i^^ < 3 X 10-^Vb, (17) 

where /xb is the Bohr magneton. The magnetic field inside the source can be so large that 
even a transition moment as tiny as lO'^^^^ can lead to a sizeable production of tau neutrino 
and anti- neutrino at the source |10]. In fact according to the models [H], it is possible to 
have 

B \ ( L \ 
109 Gauss J \10^ cm J \10-^^ fie, 
where B is the magnetic field and L is the linear size of the volume in which the magnetic 
field is as large as B. If condition (fT8|) is fulfilled, the flavor composition at source will 



> 1 
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considerably deviate from the standard 1:2:0. Let us normalize the initial flavor ratio 



such that w° = 1: 



In this section, we discuss whether by merely studying R and without theoretical prejudice, 
one can extract and w^. Actually, the possibility of deriving information on and 
from the measurement of the flavor ratio at the detector has been discussed in the literature 
|42j . Like the case of measurement of neutrino parameters, the flavor identiflcation power 
of the neutrino telescopes has not been realistically treated in the previous studies. Here 
we investigate this possibility considering realistic flavor identiflcation power of neutrino 
telescopes {i.e., studying ratio R) and taking into account various sources of uncertainty for 
the flrst time. 

There is a subtle point here. If within the lifetime of the ICECUBE (or a more advanced 
neutrino telescope) neutrino flux from a GRB in a close-by galaxy (a ~ 5 Mpc far away 
galaxy) is detected, the statistics will be high enough to extract information on the flavor 
composition of the flux from this individual source. For sources located at cosmological 
distances (> 100 Mpc), the best can be done is to combine the information from different 
sources. Each source may be different and emit neutrino flux with a different flavor compo- 
sition. From the "average" i?, only "average" values of and over these sources can be 
derived. 

Taking into account the relevant uncertainties in the input parameters, we look for values 
of i(7g : w° : that are consistent with R = R± AR. To perform this analysis, we take 
6*13 = 0. For any other value of ^13, the same analysis can be repeated. The results are robust 
against the variation of ^13 within the present bound. If ^13 = 0, by the time that enough 
cosmic neutrinos are collected, the Daya-Bay |13] and Double-Chooz [35] experiments can 
set the bound sin^ ^13 < 0.003. We vary sin^ ^13 between zero and 0.003. In this case, there 
is no hope of measuring 6 so we allow 6 to vary between and 27i. We take the energy 
spectrum to be of form E'"^ and assume that its power-law behavior will be established and 
the spectral index will be measured with 10 % precision. Again, we vary M^^/My^ within 
[0,1]. 

In Fig. ([7]), we consider two possibilities: (i) the standard case with : w° : w° = 0.5 : 
1 : leading to ^ = 2.53 (see Fig. H-a); (ii) the case of stopped muons with : w^^ : = 
0:1:0 yielding R = 3.20 (see Fig. H-b). From these flgures we observe that with a precision 
of AR/R = 7%, these two scenarios can be easily discriminated. These two can also be 
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Figure 7: Points in the {w'^,w^) plane consistent with R = R i: AR. The ratios are normahzed 



such that w?. 



1. The true values of (w^jW^j.) are denoted by Points displayed by dots and 
triangles respectively correspond to AR/R = 7% and AR/ R = 1%. In drawing this figure we have 
varied sin^ 6*13 G (0,0.003), 5 E (0,27r), a G (1.8,2.2), and TVp^TV^^ G (0,1). Drawing Fig. (a), we 
have taken R = 2.53 which corresponds to the standard picture with = 1/2 and = 0. In case 



Fig. (b), we have set R = 3.2 which corresponds to the stopped muon scenario with w\ 
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discriminated from the scenario in which the neutrino production mechanism is n ^ peue 
{i.e., : 1(7° : w° = 1 : : 0). When we restrict the analysis to = {i.e., the case 
without exotic neutrino properties) from these figures we observe that the measurement of 
R stringently constrains : which in turn sheds light on the production mechanism. 
However, once the assumption of is relaxed, a wide range of : w^^ : can be a 
solution. For example, the exotic case of : w° : w° = : : 1 leads to the same value of 
R as the stopped muon scenario. 

The input for Fig. [H] is the same as that for Fig. [7] except that in Fig. [8] Wg : w° : 
= A : 1 : 1 . For this flavor ratio, the central value of is ^ = 1.96 . Notice that with 
AR/R, this exotic flavor ratio can be discriminated from the two standard cases that we 
have mentioned. That is Wg : w° : w° = 0.5 : 1 : or Wg : w° : w° = : 1 : are not 
solutions for = 1.96(1 ± 7%) . 

7 Conclusions and Discussions 

Under the assumption that the initial flavor ratios at the source were : w^^ : = 1 : 2 : 0, 
we have studied the possibility of deriving ^13 and/or 6 from cosmic neutrinos taking into 
account various uncertainties. 

ICECUBE and other neutrino telescopes that are going to collect neutrino events in the 
energy range 100 GeV < E„ < 100 TeV will be sensitive only to two types of events; i.e., 
shower-like and yU-track events. The ratio of these two, R, yields only one piece of information 
on the mixing parameters. Under the assumption of CP conservation, cos 5 = ±1, we have 
discussed the possibility of extracting S13 from the measurement of R. We have found that 
for cos 5 = 1, the sensitivity of R to s^g is very mild. For cos 5 = 1, the derivation of S13 
from R would require measurement of R with a precision better than 2 % which does not 
seem achievable. In the case of cos 5 = —1, as 5^3 varies between zero and the present upper 
bound, R changes by 10 % which is in principle resolvable by ICECUBE [20] . 

We have found that a 10 % uncertainty in the energy spectrum (to be precise, 10 % 
uncertainty in the spectral index in Eq. ( fT3|) ) is a major source of error in the derivation 
of S13. However, for s^g > 0.02, the solutions with cos 5 = 1 and cos 5 = —1 can be 
distinguished despite a 10 % uncertainty in the spectral index. We have also studied the 
effect of a deviation from the power-law spectrum. Our conclusion is that in order to derive 
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Si3 from i?, a precision better than 5 % in the measurement of the energy spectrum is 
required. 

We have also studied the effects of the variation of the neutrino to anti-neutrino ratio 
on R. We have found that when My^/My^ (see Eq. f fT3|) for the definition) varies between 
and 1, R changes by 5 % which is comparable to the effect of S13. Unfortunately with the 
current techniques, it is not possible to measure Ny^/Nu^ in the energy range 100 GeV < 
Ey < 100 TeV so this source of uncertainty cannot be eliminated by measurement and one 
should rely on the models to predict the value of Nu^/Ny^. 

The uncertainty in the neutrino nucleon cross section, cTi^at, induces relatively large uncer- 
tainty in the evaluation of the shower-like and /i-track events. However, when we take their 
ratio, the uncertainties cancel each other out so the derivation of the neutrino parameters 
from R does not suffer from the uncertainty in (Tpn- 

We have then studied the possibility of deriving 5 from the cosmic neutrino flavor compo- 
sition assuming that S13 will be measured by other experiments with a reasonable precision. 
Since the dependence of the oscillation probabilities of the cosmic neutrinos on 5 is through 
cos 5, there is a symmetry under 5 ^ 27r — 5. On the other hand, the sensitivity of the 
CP-odd combination Pi^u^ Ue) — P{^n — > ^e) (which is proposed to be measured by neu- 
trino factory or superbeam facilities) to 5 is through sin 5 and would therefore suffer from 
a degeneracy under 5 ^ -k — 5. To resolve the latter degeneracy, it is suggested to employ 
various baselines and/or study the energy dependence of the oscillation probability [45j. 
Derivation of cos 5 from the cosmic neutrino flavor composition can be considered as an 
alternative method to solve this degeneracy. We have studied the effects induced by the 
error in the mixing angles (^12, ^23 and 6*13) and the measurement of R as well as by the 
uncertainties in neutrino-antineutrino ratio and the energy spectra. We have found that 
even with a precision of 1% in i?, CP cannot be established. That is even for the maximal 
CP- violation {5 = vr/2), 6 = cannot be ruled out. However, in this case, 5 = vr can be 
excluded provided that AR/R is reduced to less than 1%. 

Within the SM of the particles, various mechanisms can deviate w'^/w^ from 2 but within 
the energy range 100 GeV < < 100 TeV, still remains much smaller than and 
w^. The conditions in the source of cosmic neutrinos are so extreme that the beyond SM 
properties of neutrinos can play a role to significantly distort the initial flavor ratio of the 
cosmic neutrinos. In particular, a nonzero fire or firfi transition moment close to the present 
bound can lead to w° ~ w° ~ w^. In the second part of the paper, we have relaxed the 
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assumption : : = 1 : 2 : and have studied the possibihty of extracting the 
initial ifg : : w'^ from the cosmic neutrino data. We have found that with a precision of 
AR/R = 7%, one can discriminate between the standard picture with : w° : w° = 1 : 2 : 
and the stopped muon scenario with : w° : w° = : 1 : 0. When based on theoretical 
prejudice, we restrict the analysis to w° = 0, we find that w^/w'^ can be constrained with 
reasonable accuracy but relaxing = will open up the possibility of different solutions. 

We have also enumerated a number of other effects that can be comparable to that of 
nonzero S13 but have been overlooked in the literature. Calculating these effects requires the 
knowledge of details of the detector and the shape of the neutrino spectrum. Uncertainty 
in this knowledge will lead further uncertainty in the determination of 6 and ^13. As we 
demonstrated in the present paper, even in the absence of these effects, the uncertainties are 
too large to allow for the determination of 6 and 6*13. Potential uncertainties in these effects 
will further confirm this conclusion. For the purpose of establishing substantial deviation 
of : w° : from 1:2:0, these effects have to be taken into account however the 
uncertainties in the evaluation of these effects are not expected to be so large to change 
our positive conclusion in the second part of this paper. Throughout this paper, we have 
assumed that the propagation of the cosmic neutrinos from the source to the detector is 
governed by the standard oscillation formula. The effects of a deviation from the standard 
oscillation probability will be presented elsewhere. 
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A Calculation of f[Er^E^) 

In this section we calculate the function f{ET-,E^) which is the probability density of the 
emission of a muon with energy in the decay of a r lepton (r /iz/^z/,-) with energy E^ 
(see Eq. ([8]) for the definition). In the rest frame of r, the partial decay rate of an unpolarized 
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r is given by the following well-known formula (see [16 



' '^-'S] (19) 



V dE'^dVL' nml \ 3m^ 



where the effects of m^/m^ <^ 1 are neglected. Quantities in the rest frame of the decaying 
r lepton are denoted by a prime. In the rest frame of the r lepton, < -E^ < The 
number of emitted muons in certain direction within the solid angle dQ' and with energy in 
the interval [E'^, E'^ + dE'^ is a Lorentz invariant quantity: 

1 d^T 1 d'^V 

dE^dQ = - " ' dE'dQ'. (20) 



TdE^dVL^ ^ T'dE'dVL' 



From this equality, we obtain 



1 d'T _ _ 12 



dE„d^ 



TdE„dn ^ vrm! 



4 



7(1 -/?cos^)E 



-f{l- Pcos9)EldE^smeded(l), (21) 



where 7 = Er/nir and /5 = a/1 — 1/7^- The ^-axis is taken along the direction of motion of 
r. The quantities E^, 6 and take values in the following intervals 

o<0<27r, o<£;^< ^(i + /5), Q<e<e^a. (22) 



where 

, 1 / m. 
c'max = arccos max < — I 1 



13 \ 2^E^ 
By integrating over 6 and in Eq. fl2T]) . we obtain 

27r I'Omax (^^r 



(23) 



/(E.,E,)^-^-§ + i|^. (25) 



In the limit (3 1 (or equivalently, 7^1 

^Er E^ ' 3^4' 
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